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>■ I A new population of extended, luminous globular clusters has recently been 

discovered in the outskirts of M31. These objects have luminosities typical of clas- 
sical globular clusters, but much larger half-light radii. We report the first results 
from deep ACS imaging of four such clusters, one of which is a newly-discovered 
_ example lying at a projected distance of ~ 60 kpc from M31. Our F606W, 

I F814W colour-magnitude diagrams extend ~ 3 magnitudes below the horizontal 

^ I branch level, and clearly demonstrate, for the first time, that all four clusters 

are composed of > 10 Gyr old, metal-poor stellar populations. No evidence for 
^ I multiple populations is observed. From a comparison with Galactic globular clus- 

O I ter fiducials we estimate metallicities in the range —2.2 < [Fe/H] < —1.8. The 

observed horizontal branch morphologies show a clear second parameter effect 
between the clusters. Preliminary radial luminosity profiles suggest integrated 
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^ ■ ular cluster luminosity function. Our results confirm that these four objects are 

bona fide old, metal-poor globular clusters, albeit with combined structures and 



2 ' magnitudes in the range —6.6 < My < —7.7, near the median value of the glob- 
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luminosities unlike those observed for any other globular clusters in the Local 
Group or beyond. 

Subject headings: galaxies: individual (M31) — galaxies: star clusters — globular 
clusters: general 



Introduction 



As the nearest large spiral galaxy to our own, M31 is a prime target for stellar pop- 
ulations studies aimed at understanding the formation and evolution of such galaxies. An 
important subset of this work concerns the globular cluster system of M31, since the objects 
belonging to this system allow us to trace the star formation history, chemical evolution, 
and mass distribution of M31. Recently, as part of a major survey of the M31 halo with the 
Isaac Newton Telescope (INT), a search for previously unknown globular clusters, at large 
projected ra dii from the cen t er of M 31, has been conducted. Details of the INT survey may 



be found in [Ferguson et al.l (120021 ) (see also the references therein) ; spatial coverage now 
extends to more than 40 deg^. The globular cluster search has resulted in the discovery of 

th projected radii {Rp) extending to ~ 100 kpc (see 



Huxor et al. 


2004; 


Huxor 


2006) 



e.g. 

these new clusters has recently been carried out - first results from this work are described 
elsewhere ( iMackey et al.ll2006l . hereafter Paper II). 



One surprising result to arise from the globular cluster survey was the discovery of 
a previously unkno wn population of extended, luminous objects in the outskirts of M31 
( iHuxor et al.ll2005l ). Specifically, three such clusters were discovered, with 15 < i?p < 35 
kpc. These clusters are notable because they have large half-light radii (r/i), but unlike 
known globular clusters with similar structures, also have luminosities near the peak of the 
globular cluster luminosity function (GCLF). As such, they are dissimilar to any other known 
clusters in the Milky Way, M31, or anywhere else, and begin to fill in t he gap in paramete r 
space between classical globular clusters and dwarf spheroidal galaxies (IHuxor et al.ll2005l ). 



Subsequent to the discovery of these three clusters, a further survey of the southern 
quadrant of the outer M31 halo, and down the minor axis of M31 to M33, has been cond ucted 
with MegaCam on the Canada- France-Hawaii Telescope (see e.g., iMartin et al.ll2006l ). An 
additional extended cluster was located wi th these obs e rvatio ns, at a projected radius of 
~ 60 kpc. Following the nomenclature of iHuxor et al.l (120051 ). we label this object EC4. 
Its coordinates are listed in Table [1] This is one of the most remote M31 globular clusters 
currently known, after two of the newly-discovered classical globular clusters at Rp ~ 78 and 
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100 kpc (see Paper II), and the new cluster of iMartin et al.l (120061 ) . at Rp ~ 116 kpc. 

In this Letter, we report on the first resuhs from imaging of the four unique luminous, 
extended M31 clusters with the Advanced Camera for Surveys (ACS) on-board HST. 



2. Observations and data reduction 

Our observations were obtained with the ACS Wide Field Channel (WFC) under HST 
program GO 10394 (P.I. Tanvir), between 2005 July 10 - 2005 July 23. Each cluster was 
imaged three times in the F606W filter and four times in the F814W filter, with small dithers 
between subsequent images. Typical total integration times were 1800s in F606W and 3000s 
in F814W. Each cluster was placed at the centre of either chip 1 or chip 2 on the WFC, 
to ensure that the inter-chip gap did not impact the observations. In Fig. [1] we display 
drizzled F606W images of the four clusters, confirming the nature of these objects. Each 
is clearly very extended and of low surface brightness, quite different in structure to the 
cl assical globular c lusters which have previously been observed in M31 (see e.g., the images 



m 



Rich et al.ll2005f ) 



We used the dolphot photometry software (lDolphinll2000l ). specifically the ACS mod- 
ule, to photometer our images, dolphot performs point-spread function (PSF) fitting using 
PSFs especially tailored to the ACS camera. For a given cluster, photometry was done si- 
multaneously on all the flat-fielded images from the STScI archive (both filters), relative 
to some deep reference image - we used the drizzled combination of our F814W images. 
DOLPHOT accounts for the hot-pixel and cosmic-ray masking information provided with 
each flat-flelded image, fits the sky locally around each detected source, and provides mag- 
nitudes corrected for charge- transfer efficiency (CTE) effects on the calibrated VEGAMAG 



scale of ISirianni et al.l (120051 ). A variety of quality information is listed with each detected 
object, including the object type (stellar, extended, etc), of the PSF fit, sharpness and 
roundness of the object, and a "crowding" parameter which describes how much brighter 
an object would have been had neighbouring objects not been fit simultaneously. We used 
some of this information to select only stars with high quality measurements. Specifically, 
we selected objects of stellar type, with valid photometry on all seven input images, global 
sharpness parameter between —0.3 and 0.3 in each filter, and crowding parameter less than 
0.25 in each filter. 

We isolated each cluster from the surrounding field population by imposing a radial cut. 
Two (ECl, EC3) are located, in projection, fairly close to the main body of M31, and hence 
have heavy background fields. For these objects, we set the selection radii to be small (15" 
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and 8", respectively). We found these sufficient to obtain clean cluster colour-magnitude 
diagrams (CMDs) with plenty of stars still available to define the primary sequences. For 
the two clusters at large projected radii (EC2, EC4), we used limiting radii closer to the 
expected cluster tidal radii, at ~ 25 — 30". 



3. Analysis 

Fig. [2] shows CMDs for the four clusters. Our photometry reaches ~ 3 mag below the 
level of the horizontal branch (HB), to a limiting magnitude of mpeoew ~ 28. Each cluster 
exhibits a narrow, steep red-giant branch (RGB), and clearly delineated horizontal branch. 
Three of the clusters possess HBs extending to the blue, including broadened regions with 
colours in the range 0.2 < mFeoew — '"^F814W ^ 0.5 which are suggestive of RR Lyrae stars 
imaged at random points during their cycle. Taken together, these features are indicative of 
old (> 10 Gyr), metal-poor stellar populations. 

We observe no sign of young or intermediate-age populations in these clusters, such as a 
blue plume or luminous AGB stars above the RGB tip. Comparing CMDs for the inner and 
outer parts of each cluster does not re veal any radially dep endent differences such as those 



commonly seen in dSph galaxies (e.g.. iHarbeck et al.ll200ll ). Furthermore, the widths of the 



four RGBs above the HB level are comparable to each other and to those for our classical 
globular clusters (see Paper II), and are consistent with the photometric errors. Hence, we 
see no evidence for multiple stellar populations within our four clusters. This argues against 
merger formation scenarios except those in which the parent clusters were of similar age and 
metallicity, as well as scenarios involving stripped dSphs. 

In order to obtain photometric metallicity estimates from our CMDs, we us ed the cluster 



f iducia l sequences in the ACS/WFC F606W and F814W bands presented by iBrown et al. 



( 120051 ). Because our four targets are distributed across a large area of sky and a wide 
variety of projected radii from M31, it is prudent not to assume a single distance modulus 
or reddening value for all. Rather, we adopted a procedure to simultaneously constrain 
the metallicity, distance modulus, and reddening for each cluster. We chose to register the 
fiducials to our cluster CMDs using the F606W magnitude of the HB and the colour of the 
red-giant branch (RGB) at the HB level. If the metallicity of the reference fiducial matches 
that for a given cluster, then when registered using these two points, the fiducial will closely 
trace the upper RGB of the cluster. If the reference fiducial is of incorrect metallicity, the 
fiducial and CMD will deviate on the upper RGB. Bounding the cluster RGB with more 
metal-rich and metal-poor fiducials allows the true cluster metallicity to be estimated. 
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The absorption in ea ch ACS bandpass i s a non-hnear function of both effective tem- 
perature and E{B — V). iBrown et aLl (120051 ) therefore provide two coefficients {a, (3) per 



bandpass at each point on their fiducial sequences, allowing the absorption to be calculated 
given E{B — V). For example, in F606W the absorption Apsoew is: 



^F606W = OfpmswE{B — V)+ (3Ymsw[E{B — V)Y 



and similarly in F814W. Following iMcConnachie et al.l (120051 ). M31 has distance modulus 
(m — M)o = 24.47 ± 0.07; however, due to the line-of-sight depth of the system, the cluster 
distance moduli must be scattered about this value. We decided to test over a range of ±0.5 
mag, representing a system depth of ~ 360 kpc. For a given (m — M)o, we used Eq. [T]to solve 
for the E{B — V) which aligned the HB levels of the observed cluster and the fiducial under 
consideration. With this constrained, we calculated the difference between the colour of the 
cluster RGB at HB level, and that of the transformed fiducial; the best solution for (m — M)o 
and E{B — V) was the combination which minimised this offset. We derived uncertainties in 
these quantities by randomly selecting new HB levels for the fiducial and cluster and a new 
RGB colour for the cluster, from Gaussian distributions about our measured values (typical 
distribution widths were ±0.1 mag in the HB levels, and ±0.01 mag in the RGB colour), and 
solving for new (m — M)o and E{B — V). Iterating ten thousand times built up distributions 
from which we derived standard random errors. These are generally 0.14 mag in (m — M)o 
and 0.01 mag in E{B — V). We note that there are possibly additional systematic errors 
present for our measured E{B — V) values, due to differences in age between the M31 clusters 
and the fiducial clusters, and the fact that the colour of the RGB at the HB level is mildly 
age-sensitive. At most, these errors would amount to ~ 0.02 mag more than quoted in Table 
[H corresponding to the interval 10 — 15 Gyr. 

With the fiducials aligned, it was straightforward to infer the most likely value of [Fe/H] 
(and hence the best values for (m — M)o and E{B — V)) for a given cluster. The fiducial 
alignments are presented in Fig. [3l and the numerical results in Table [H Typical errors in 
[Fe/H] are approximately ±0.15 dex. From Fig. [3l it is clear that all four clusters are metal- 
poor objects. Three have an RGB locus intermediate between those for NGC 6341 ('M92) and 
NGC 6752, while the fourth (EC3) matches closely that for NGC 6341. iBrown et al.l J2005h 
adopt [Fe/H] = —2.14 and —1.54 for these two clusters, respectively; the listed metallicities 
in Table [1] reflect these values. 



For comparison, we consider values of E{B—V) from the extinction maps of lSchlegel et al 



(119981 ) (see Tabled]). The agreement is close, indicating the vali dity of our method. W e note 
a tendency to over-estimate E{B — V) by 0.01 — 0.02 mag; iBarmby et al.l (l2000l ) found 
a similar effec t in th eir large study of M31 globular clusters. Adopting E{B — V) from 
Schlegel et al.l (119981 ) tends to increase (m — M)o by a few hundredths of a magnitude and 
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does not noticeably affect tlie derived [Fe/H]; liowever tlie alignment for EC4 is noticeably 
inferior. Our derived distanc e moduli for the tAvo clust ers with smallest Rp are ~ 24.45, 
matching the M31 value from iMcConnachie et al.l (120051 ). EC4 may lie somewhat closer to 
us than the bulk of M31; however this is unsurprising given its large Rp. 

From our CMDs and estimated [Fe/H], it is clear that there is a strong second parameter 
effect among this set of clusters. Deriving accurate HB-morphology parameters requires 
variability information, completeness test s, and a treatment of photometric blends at the red 
end of the HB (see e.g., iRich et al.ll2005l ). and is therefore beyond the scope of the present 
letter. We defer such calculations to a future paper, currently in preparation. Even so, 
examining the three clusters with very similar [Fe/H] ~ —1.84, the second parameter effect 
is easily discernible by eye. EC2 has an exclusively red HB, while EC4 has an intermediate 
HB and ECl has a mostly blue HB. The CMD for EC2 is strikingly similar to that for the 
SMC globular NGC 121, which has comparable [Fe/H]. This cluster is ~ 2 — 3 Gyr younger 



than the oldest Galactic and LMC globular clusters (see e.g.. lShara et al.lll998l ). suggesting 
that a similar age spread may be present among our M31 clusters. 

Finally, we obtained estimates of structural parameters and total luminosities by means 
of integrated (aperture) photometry. This technique ensures unresolved cluster light is 
counted, and obviates any present need for lengthy completeness tests. Using the CMDs 
we estimated the largest apertures (rmax) within which field star contamination does not 
add significantly to the integrated cluster luminosity. For a given cluster we used this aper- 
ture to obtain the total luminosity in F606W and F814W (background galaxies and bright 
field stars were masked) , then used these to transform the F606W value into Johnson V 
using the prescription in ISirianni et al.l (120051 ) . We then calculated My using our measured 
E{B — V) and (m — M)q. These luminosities are lower limits, since r^ax is almost certainly 
smaller than the tidal radius in all clusters. For EC2 and EC4 the amount of missing light 
is only ~ 0.1 mag, since rmax ~ 75 pc. However, for ECl and EC3, which have small r^ax, 
the shortfall is significant. Based on profiles for EC2 and EC4, the true luminosity of ECl 
may be ~ 0.3 mag brighter than listed, and that for EC3 ~ 0.7 mag brighter. 

With My calculated, we sampled the luminosity profile of each cluster with a variety 
of apertures to determine the half-light radius r^. Subtracting the luminosities within con- 
secutive apertures allowed preliminary surface brightness profiles to be constructed. From 
each, the core radius Tc (at which the surface brightness is half its central value) was esti- 
mated. Table [1] lists Vc and Vh, converted to parsecs using our measured (m — M)q. Typical 
uncertainties are ±15% in each. 



Taken together, the above results confirm the original conclusion of iHuxor et al.l (120051 ) 
that these four objects strongly resemble classical globular clusters in terms of their old. 
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metal-poor stellar populations, but are quite different structurally. Specifically, their lumi- 
nosities sit around the peak of the GCLF, and as such they are much bright er than similarly 
exten ded objects found in the Milky Way (including the n ew cluster of iBelokurov et al. 



20061 ). the LMC, or externally (such as the "faint fuzzies" of iBrodie &: Larsenll2002l ). They 
thus encroach on the gap in parameter space between classical globular clusters (which are 
dark matter free) and dwarf spheroidal galaxies (which are dark matter dominated). Fur- 
ther study of this unique population will shed light on the transition region between globular 
clusters and dwarf spheroidals, and provide important clues to the formation and evolution 
of these clusters, which may be strongly linked to that of M31. 
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Table 1. Observed properties of four luminous, extended globular clusters in the outskirts of M31 



Identifier'* 


RA 


Dec 


Rp 


(m - M)o 


E{B - V) 


E{B - V) 


[Fe/H] 


rc 






^max 




(J2000.0) 


(J2000.0) 


(kpc) 




(meas.) 


(lit.) 




(pc) 


(pc) 




(") 


ECl 


00'>38™19!5 


+41°47'15'.'0 


13.3 


24.42 ±0.15 


0.08 ±0.01 


0.07 


-1.84 


20.4 


35.4 


-7.4 


16 


EC2 


00''42"55'!1 


+43°57'28'.'5 


36.8 


24.39 ±0.14 


0.10 ±0.01 


0.09 


-1.84 


13.2 


29.5 


-7.0 


20 


EC3 


00''38'°04!5 


+40°44'37'.'5 


14.0 


24.47 ±0.14 


0.07 ±0.01 


0.07 


-2.14 


15.1 


32.3 


-7.0 


9 


EC4 


00'*58"15!5 


+38°03'01'.'l 


60.0 


24.31 ±0.14 


0.08 ±0.01 


0.05 


-1.84 


23.8 


33.7 


-6.6 


20 



^Nomenclature follows that of Huxor et al. ( 20051 ). with the addition of a new cluster, EC4. 

''Integrated luminosities are lower hmits, as discussed in the text. In particular, ECl and EC3 may be ~ 0.3 and ^ 0.7 mag brighter 
than listed, respectively. 
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Fig. 1. — Drizzled ACS/WFC F606W images of the four extended luminous M31 globular 
clusters. Each thumbnail has dimensions of 25" x 25". 
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Fig. 2. — CMDs for the four extended, luminous M31 globular clusters. Photometry has 
been selected from the full ACS/WFC fields by imposing radial limits of, respectively, 15", 
25", 8", and 30" from the cluster centers. Typical photometric errors from dolphot are 
indicated. 
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Fig. 3. — Results of fitting the Galactic globular cluster fiducials from iBrown et al.l ( l2005l ) 
to the observed CMDs. All four clusters are metal-poor, and hence have RGBs bracketed 
by that for NGC 6341 (M92) at [Fe/H] = -2.14 (sohd lines), and that for NGC 6752 at 
[Fe/H] = —1.54 (dashed lines). EC3 is the most metal poor cluster, with the M92 fiducial 
providing a close fit. The horizontal dotted lines represent the adopted HB levels. The 
numbers in the upper left indicate the best fitting (m — M)o and E{B — V) for the rnarked 
fiducials. Photometry has been transformed to STMAG to match the iBrown et al.l (120051 ) 
fiducials. 



